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We present generalized scattering theory for multi-terminal spin transport in systems with broken 
SU(2) symmetry either due to spin-orbit interaction, magnetic impurities or magnetic leads. We 
derive equation for spin current consistent with charge conservation. It is shown that resulting spin 
current equations can not be expressed as difference of potential pointing to non applicability of 
linear response for spin currents and as a consequence equilibrium spin currents(ESC) in the leads 
are non zero. We illustrate the theory by calculating ESC in two terminal normal system in presence 
of Rashba spin orbit coupling and show that it leads to spin rectification consistent with the non 
linear nature of spin transport. 
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Spin transport has emerged as an important subfield 
of research in bulk condensed matter system as well in 
mesoscopic and nano system 1]. In macroscopic systems, 
the very definition of spin currents is still debated due to 
non conservation of spin in presence of SO interaction^ 
. On the other hand in mesoscopic hybrid system since 
current is defined in the leads where SO coupling is 
absent, therefore, it has been assumed that Landauer- 
Biittiker formula for charge current [1] (Eq. ([9]) in this 
manuscript), which determines current in leads in terms 
of applied voltage difference multiplied by total transmis- 
sion probability, can be straight away generalized for spin 
currents by replacing total transmission probability with 
some particular combination of spin resolved transmis- 
sion probabilities This simple generalization has 
been widely used in the literature to study spin depen- 
dent phenomena in nanosystems0, [^, @. The Landauer- 
Biittiker formula in its widely used form has inbuilt cur- 
rent conservation ,z.e., total current is divergence less 
(div j=0) which follows from basic Maxwell equations of 
electrodynamics. Physically this implies that total cur- 
rent has neither sources nor sinks this would be true for 
spin currents as well if spin is conserved. However, in 
presence of spin-orbit interaction, magnetic impurities or 
non-coUinear magnetization in leads, spin is no longer a 
conserved quantity, hence the spin currents can not be 
divergence less. Therfore a straight forward generaliza- 
tion of Landauer-Biittiker charge current formula to spin 
currents can not be correct for spin non-conserving sys- 
tems. 

In view of the above discussion in this work we de- 
velop a consistent scattering theoretic formulation of cou- 
pled spin and charge transport in multiterminal systems 
with broken SU(2) symmetry in spin space following 
Biittiker's work on charge transport [3i|. The SU(2) sym- 
metry in spin space can be broken due to either SO inter- 
action, magnetic impurities or non-coUinear magnetiza- 
tion in leads . Our analysis provides a correct gener- 
alization of Landauer-Biittiker theory for spin transport. 
In particular we derive a spin currents equation (Eq. ([5]) 



in this manuscript) consistent with charge current conser- 
vation. However, the resulting spin current equation can 
not be cast in terms of spin resolved transmission and re- 
flection probabilities multiplied by voltage difference as is 
the case for charge current(Eq. @ in this manuscript)^. 
Therefore, equilibrium spin currents are generically non 
zero. Moreover, it implies that linear response theory 
with respect to electric field is not applicable to the spin 
currents (equilibrium as well non- equilibrium). Thus spin 
currents are intrinsically non-linear in electrical circuits. 
However, this is not surprising since linear response is 
valid for thermodynamically conjugate variable. In an 
electrical circuit thermodynamically conjugate variable 
to electric field is charge current not the spin currents. 
We illustrate the theory by calculating ESC analytically 
for two-dimensional electron system with Rashba SO in- 
teraction in contact with two unpolarized metallic con- 
tacts. Our analytical formula for ESC clearly demon- 
strates that it is transfer of angular momentum per unit 
time from SO coupled sample to leads where SO inter- 
action is zero. Therefore, it is truly a transport current 
in contrast to equilibrium spin currents in macroscopic 
Rashba medium(see ref.jl, Q). 

To formulate scattering theory for spin transport we 
consider a mesoscopic conductor with broken SU(2) sym- 
metry in spin space connected to a number of ideal mag- 
netic and non-magnetic leads (without SO interaction) 
which in turn are connected to electron reservoirs. To 
include the effect of broken SU(2) symmetry, it is neces- 
sary to write spin scattering state in each lead along local 
spin quantization axis. For magnetic leads local magneti- 
zation direction provides a natural spin quantization axis 
which we denote in a particular lead a by rha{'daT'^a) 
where "da and (pa is polar and azimuthal angle respec- 
tively. For nonmagnetic leads since there is no preferred 
spin quantization axis, hence we choose an arbitrary spin 
quantization axis u{9, (p) which is same for all nonmag- 
netic leads. Thus the most general spin scattering state 
in lead a which can be either magnetic or nonmagnetic 
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is given by, 



^':^(r,t)= / dE y , 

«„(i?)e^^^"(^)- + fo^„(i?)e-''=="(^)-) (1) 

where ^an{'>'±) is transverse wavefunction of channel n 
and Xa(''') is corresponding spin wave function along 
chosen spin quantization axis, u or nia such that S • 
iix((7) = (crS/2)x(o-) or S ■rna'^{<j)—{ah/2)ip{a) with <t = 
(7 ( (T=±l, representing local up or down spin compo- 
nents) for nonmagnetic and magnetic leads respectively. 
Here S = {h/2)a is a vector of Pauli spin matrices and 
is number of channels with spin a in lead a. The rela- 
tion between spin dependent wavevector k^^i^) ^^'^ en- 
ergy E is specified by, E = [fi2fc^„^/2m -I- £„„ -I- a A^] , 
where is energy due to transverse motion, is 
stoner exchange splitting in the magnetic lead a. The 
stoner exchange splitting is zero for nonmagnetic leads. 
The operators aj^„ and 6J^„ are annihilation operator for 
incoming and outgoing spin channels in lead a and are 
related via the scattering matrix. 
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The scattering matrix elements S^'^.p^ provides scatter- 
ing amplitude between spin channel na' in lead /3 to spin 
channel ma in lead a. These scattering matrix elements 
will be function of energy E as well angles, d and ip. 
The angular dependence of scattering matrix elements 
on polar and azimuthal angle arise due to broken SU(2) 
symmetry. Note that for noncoUinear magnetization in 
leads and in absence of SO interaction and magnetic im- 
purities, the angular dependence is purely of geometric 
origin and is related to the angular variation of various 
magnetoresistance phenomena |10| . 

The current in spin channel a along longitudinal direc- 
tion i (through a cross section of lead a) and the local 
spin quantization axis u is defined as. 
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can take on values icr corresponding to two spin pro- 
jections along local spin quantization axis. The second 
term of Eq.(j4]) can be written explicitly in terms of spin 
resolved reflection and transmission probabilities as, 

En^a a Qoa' \ ^ a\ a' a naa' 

13a' -.nin 
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Where i?^^ and T^p are spin resolved reflection and 
transmission probability in the same probe and between 
different probes respectively. In Eql5]on right hand side 
spin resolved reflection and transmission probabilities are 
summed over all possible input modes for a fixed output 
spin mode a in lead a. Because partial scattering matrix 
in spin subspace is not unitary due to non conservation of 
spin hence this summation need not to be equal to num- 
ber of spin a channels in lead a, i.e. N^, rather it can 
have any value lying between zero and N^. To determine 
in terms of spin resolved reflection and transmission 
probabilities, consider a situation where current is in- 
jected from reservoir only in spin channels a in lead a. 
In this case charge conservation requires that this current 
should leave the spin channel a through all other possible 
channels in the same lead as well in differing leads, which 
implies, 
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Substituting for from Eq. ^ into Eq. ([H we get an ex- 
pression for spin current in terms of creation and annihi- 
lation operators. On the resulting expression we perform 
quantum statistical averages and after a lengthy algebra 
we obtain following expression for average current in spin 
channel a (for brevity of notation we suppress the super- 
script X written in Eq. 
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Where fp = l/exp[{E — ^p)/kT] + 1 is Fermi distribu- 
tion function with chemical potential /i/j and the pre- 
factor g equals ah/2. The summation over a' in Eq. Q 



As we can see that Eq.® differs from Eq.([5]) in a sub- 
tle way and are not equal because in general spin re- 
solved transmission or reflection probabilities can not be 
related among themselves by interchanging spin indices, 
i.e., T;^'f ^ T^^ and i?^""^ ^ R'^" ( we wifl discuss 
constraints due to time reversal symmetry below). If we 
demand that sum in Eq.® also equals to then it 
would imply spin conservation which is incorrect in pres- 
"^'nce of spin flip scattering or broken SU(2) symmetry. 
The inadvertent use of this charge conservation sum rule 
for spin degrees of freedom in Ref. |Tl| has led to 

incorrect spin current equation. Though the partial scat- 
tering matrix in spin subspace is not unitary, however, 
the full scattering matrix is unitary,i.e., SS'^ = S'' S = I, 
therefore, if we sum over a also in Eq.® or Eq.® then 
it should give total number of channels in leads a, i.e. 
N = + N~'^ , and as a result we get the following sum 
rule for total transmission probability, 
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where Tap is total tranmission probability. 

The net spin current flowing in lead a is deflned as 
^iia — i^iia) + i^iil) while the net charge current flowing 
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is given by sum of absolute values, i.e., i-^Za) I + I 

i^iia) I "^ith pre- factor g replace by the electronic charge 
e in Eq. (g]). Using Eqs.( l^),® in Eq.® we obtain net 
spin and charge current as. 
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Equation ([5]) is the central result of this work. We stress 
that Eqs. (0) and (0) are valid under most general condi- 
tions as we have not made any assumptions about sym- 
metries of the scattering region. It is instructive to note 
that in general T;^'/ + T^'" and R"^-" ^ R'^" there- 
fore, spin current equation can not be simplified further 
and written in terms of difference of Fermi function mul- 
tiplied by transmission or reflection probabilities as is 
the case for charge current in Eq. ([9]) which is stan- 
dard Landauer-Biittiker r esult'a^. Hence the spin current 
given by Eq. ([8]) will be nonzero even when all the leads 
are at equilibrium, i.e., fa{E,^ia) — fiE,^),\/a, where /i 
is equilibrium chemical potential. For sake of complete- 
ness we mention that equilibrium charge current vanishes 
as is evident from Eq. ([9]). The preceding discussion im- 
plies that linear response for spin currents is not appli- 
cable in an electrical circuit where external perturbation 
is applied voltages which is conjugate to charge currents 
and not to the spin currents as discussed in introduction. 
Therefore, the most widely used equation for spin cur- 
rent, see Ref.fi', "g"] obtained by a generalization of charge 
current Eq. ([9]) has to regarded as incorrect. In view of 
this the theoretical study of spin dependent phenomena 
in mesoscopic systems needs to be re-investigated. 

We can gain further insight into spin current by con- 
sidering non-equilibrium situation such that the chemi- 
cal potential at the different leads differ only by a small 
amount so that we can expand the Fermi distribution 
function around equilibrium chemical potential fi as, 
f(3{E,fii3) = f{E,n)-{-{-df/dE){n0-ii). In this case we 
can immediately notice from Eq. ^ that total spin cur- 
rent in non-equilibrium situation will have equilibrium as 
well non equilibrium parts of spin current. For ESC the 
full Fermi sea of occupied levels will contribute. There- 
fore even in non-equilibrium situation ESC cannot be 
neglected. 

Equilibrium spin currents in time reversal symmetric 
two terminal system: In time reversal symmetric systems 
spin resolved transmission and reflection probabilities in 
Eq. ^ obey following relations i.e., i?^^ = i?^^ and 



Ta^ [3] ■ In this case the spin currents Eq. 



nals by L and R respectively), 
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above equation gives spin current in Left terminal. Spin 
current in right terminal are obtained from the same 
equation by interchanging L <^ R. On right hand side 
in Eq. pop . a and —a refers to up and down spin states 
along u. From the above equation and previous discus- 
sion it is evident that even in time reversal symmetric 
two terminal systems ESC are non zero. Incase SU(2) 
symmetry in spin space is preserved, the spin resolved 
transmission and reflection probabilities obey a further 
rotational symmetry in spin space, i.e, T^^ — T^^^'^ , 

— ^aa '^ ^-ii^ Spin flip Components are zero, which 
implies that spin currents are identically zero for all ter- 
minals as is evident from Eq. (jlOp . This conclusion re- 
mains valid even for systems without time reversal sym- 
metry as can be seen easily from Eq. (|S]). 

The expression in Eq. (jlOp can be cast in a more useful 
form as (the details will be provided in Ref . [T^ ) . 



^ j Tr„ [{F^GT^G'^-f 
{V^g''T^g'')}[a-u)]f{E,ij)dEdk. 



(11) 



In the above equation all symbols represents 2x2 matri- 
ces in spin space( in cr • m basis) and trace is taken over 
spin space. Where Ta^p represents broadening matrices 
due to contacts, G''^"^ are retarded and advanced Green 
function and g'"'" is a off diagonal matrix in spin space 
defined as g''-'* = [{0, G^'^^}, {G!:"^, 0}]. First term in 
Eq. (|lip corresponds to spin resolved transmission while 
the second and third term give spin resolved reflection 
probabilities as required by Eq. (fTO|) . Notice that the 
above formula can not be simply written in terms of 
transmission matrix and it is reminiscent of the charge 
current formula for interacting system derived in Rcf.[12j]. 
In our case this happens for spin current because in pres- 
ence of SO interaction spin can not be described as a 
non-interacting object. 

Equilibrium spin currents in two terminal Rashba sys- 
tem: We now apply Eq. (fTT|l to study ESC(zero tempera- 
ture) in a finite size Rashba sample of length L, contacted 
by two ideal and identical unpolarized leads. The Hamil- 
tonian for two-dimensional electron system with Rashba 
SO interaction and short-range spin independent disor- 
der is H = S2fc2/(2TO*)J -F Xso{(Jxky - (Tykx) + U{x,y), 
where Xso is Rashba SO coupling strength, U {x, y) is the 
random disorder potential and / is 2 x 2 identity ma- 
trix. Neglecting weak localization effects, the disorder 
averaged retarded Green function including the effect of 
leads is given bvp^. 



G''^''{E,k) = 



E 



2m 



iT]{k) + Xso{o- 



further simplifies to (here we denote left and right termi- 



[{E-'^+lV{k)f-{\sokY 



(12) 
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with r]{k) = {2j{k) + where T{k) is momentum 

relaxation time due to elastic scattering caused by im- 
purities and 7(fc) broadening due to leads. The contact 
broadening matrices in Eq. are diagonal in spin space 
and defined as Ti^a = [{7(fc), 0}, {0, 7(fc)}]. Physically 
significance of j{k) is that it represents h/2 times the 
rate at which an electron placed in a momentum state 
k will escape into left lead or right lead, hence as a first 

approximation wc can write, 7(fcj — — |^ = — ' , 
where (f) is angle with respect to x axis. The impurity 
scattering time can be approximated as w , where 
lei is elastic mean free path. With these inputs we can in- 
tegrate Eq. (jlip over transverse momentum (multichan- 
nel case) and energy to obtain an analytical expression 
for equilibrium spin current. We find that ESC with 
spin parallel(antiparallel) to the z or z axis and flowing 
to the X direction vanishes in both leads(/f ^ = azVx=Q, 
= <JxVx—Q)- The ESC with spin parallel(antiparallel) 
to the y axis are nonzero and given by, 



^L,+v - ^R-y 



(2 + (^)'). 



(13) 



In left lead ESC is polarized along +y direction and flows 
outwards from sample to lead, i.e., along —x direction, 
while in the right lead ESC is polarized along —y and 
flows outwards from sample to lead, i.e., along +x direc- 
tion. Physically this implies that spin angular momen- 
tum is generated in sample with SO coupling which then 
flows outwards in the regions where SO coupling is zero, 
i.e., the left and right leads. This implies a spin rectifi- 
cation effect which can only occur if the transport is non 
linear and we see that this consistent with nonlinear na- 
ture of spin currents as remarked earlier. It is important 
to note that due to ESC there is no net magnetization in 
the total system (sample-|-leads) which is consistent with 
the Kramer's degeneracy. 

We can gain a deeper understanding of the above ex- 
pression if we analyze the systems using additional sym- 
metries. The disorder averaging establishes reflection 
symmetry with respect to x(y) axis and the system has 
a symmetry related with the operator cTyRy. (a^Rx)- 
As a result the total symmetry operator (time rever- 
sal-|-reflection) for the system is Urn — It<^xRx<^yRy — 
It{io'z)RxRy, where It is time reversal operator. Under 
this symmetry operation the disorder averaged system 
is invariant and the spin current operators axVx, cTyVx 
are even while UzVx is odd. Therefore the spin current 
along z direction vanishes while in-plane spin current can 



be non zero. However as we have seen above that only 
the y component of spin current survives after integrat- 
ing over momentum. Fig.(l) illustrate the conservation 



of 



' y ^x 



under these symmetry operation. The depen- 



dence of ESC on Aso can also be inferred from symmetry 
consideration. As we can check easily that the Rashba 
SO interaction changes sign under reflection along z 
8i-Kis{Xso{crxky - aykx) ^ -XsoicTxky - aykx)). Physi- 
cally it corresponds to reversing the asymmetry of con- 
flning potential along z axis. Therefore the spin currents 



spin parallel to - - - - 
+ y direction 



flow parallel to 
-ve X direction 




flow parallel to 
' +ve X direction 



- » spin parallel to 
-ve V di rection 



FIG. 1: Fig 1; The figure illustrate the conservation of spin 
current OyV^ if th system is rotated by tt along z axis, which 
represents reflection along and y axis respectively. Under this 
transformation configuration in left and right goes over into 
each other hence the spin current remains invariant. 



(equilibrium as well non-equilibrium) can only depend on 
the odd powers of spin orbit coupling constants Aso • Ac- 
cording to Eq. ESC are proportional to A^o which is 
consistent with these symmetry consideration and ESC 
vanishes if Aso zero as expected on physical grounds. It 
is worth noting that even the local ESC in macroscopic 
Rashba medium, discussed in Ref.0 are proportional to 
Aso ^'^'^ ^ have seen this is a consequence of symmetry 
consideration. Moreover ESC are proportional to length 
of SO region because SO region acts as source of these 
currents. Note that the right hand-side of Eq. ([TH)) has 
dimension of angular momentum per unit time signifying 
that these currents are truly transport current. 

To conclude, we have derived spin current formula 
for multiterminal spin transport for system with broken 
SU(2) symmetry in spin space. We have demonstrated 
that spin currents are fundamentally different from the 
charge currents and the ESC are generically non zero. In 
view of this the spin transport phenomena in mesoscopic 
system needs a fresh look. Further it will also be inter- 
esting and desirable to study different magnetoresistance 
phenomena from the perspective of spin currents. 

I acknowledge helpful discussion with M. Biittiker and 
A. M. Jayannavar. 
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